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ABSTRACT 

We investigated the properties of the Z pattern of the low-mass X-ray binary 
GX 17-1-2 in the X-ray colour-colour diagrams (CDs) and hardness-intensity diagrams 
(HIDs), and of the power spectra as a function of position in the Z. We used all 
EXOSAT ME data on GX 17-^2, a total of 6 observations during 1983-1986. We find 
that the Z pattern in the CD does not move to within ^^2 per cent between observations 
separated by '^l day to '^1.5 years. The power spectra and the characteristics of the 
Z pattern both support the division of the Z sources into two sub-groups, GX 17+2 
being in the group together with ScoX-1 and GX 349-1-2. 

It has been noted previously that the maximum frequency of the horizontal 
branch quasi-periodic oscillations in GX 17-1-2 is lower than for CygX-2, GX5-1 and 
GX 340-1-0. We suggest that an asymmetry in the magnetic field giving rise to polar 
caps with different emission characteristics may be the origin of this. This may also 
explain the occurrence of bursts in GX 17-1-2. 

We also report the detection of two new bursts in the EXOSAT data, of GX 17+2, 
which we interpret as type I bursts. No regular pulsations are observed during these 
two bursts nor during the bursts previously detected in GX 17-1-2 with EXOSAT, with 
upper limits to the pulsed amplitude of ~1.5 per cent for pulsations with a frequency 
below 512 Hz. Three of the four bursts occurred in the normal branch, while one 
occurred in the lower part of the flaring branch. 

The Poisson level in power spectra of data processed by the on-board computer of 
EXOSAT is found to deviate from that expected from simple variable dead time effects 
in GX 17-1-2 and also in several other sources. We present corrections to the usually 
employed relations between count rates observed with EXOSAT and the Poisson level. 

Key words: accretion, accretion disks - binaries: close - stars: individual: GX 17+2 
- stars: neutron - X-rays: bursts - X-rays: stars 



1 INTRODUCTION 

Studies of the broad-band X-ray spectral and correlated tim- 
ing behaviour of the brightest persistent low-mass X-ray bi- 
naries (LMXBs) have shown that GX 17+2 is one of the six 
known Z sources (Hasinger & van der Klis 1989), which share 
many characteristics. In the X-ray colour-colour diagram 
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(CD) and hardness-intensity diagram (HID) the sources 
trace out a roughly "Z" shaped track. The limbs of the Z 
are called horizontal branch (HB), normal branch (NB) and 
flaring branch (FB), from the top limb to the bottom limb, 
respectively. The sources move smoothly through the dif- 
ferent branches, and do not jump from branch to branch. 
This one-dimensional motion through the Z is thought to 
be governed by the mass-accretion rate, M , which increases 
from the HB, through the NB, to the FB (e.g. Hasinger et 
al. 1990). The Z sources are believed to accrete mass at a 
rate near the Eddington limit (e.g. Lamb 1989). 
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In each branch, a Z source shows a characteristic power 
spectrum (Hasinger & van der Klis 1989). In the HB one 
sees intensity correlated high-frequency (15-55 Hz) quasi- 
periodic oscillations (QPO), called horizontal-branch oscil- 
lations (HBO), together with associated band-limited noise, 
called low- frequency noise (LFN), which dominates the re- 
gion from ~1 Hz to the HBO frequency. In the NB one sees 
another type of QPO at a lower frequency (4-7 Hz), called 
normal-branch oscillations (NBO). NBO have been seen to 
occur simultaneously with HBO (e.g. Hasinger et al. 1990, 
Lewin et al. 1992). In some Z sources (ScoX-1, see e.g. Di- 
eters & van der Klis [1996], and GX 17-1-2, see e.g. Pen- 
ninx et al. [1990]), the NBO merge smoothly into higher 
frequency (up to ~20Hz) QPO, called flaring branch oscil- 
lations (FBO) when the source moves onto the FB. Recently, 
kHz QPO have been discovered in the Z source ScoX-1 (van 
der Klis et al. 1996), whose properties also seem to depend 
on position in the Z. 

Two noise components have been seen to occur in all 
three branches. One is the very-low-frequency noise (VLFN) , 
dominating the region below ~lHz. VLFN generally in- 
creases towards the FB. The other component is the high- 
frequency noise (HFN), which dominates the region above 
^20 Hz. HFN generally decreases towards the FB (see Di- 
eters & van der Klis 1996, and references therein). 

It has become clear that on the basis of their spec- 
tral and timing behaviour the Z sources can be divided 
into two categories (Hasinger & van der Klis 1989, Penninx 
et al. 1991, Kuulkers et al. 1994a), the difference between 
them may be due to a difference in orbital inclination an- 
gle (Hasinger & van der Klis 1989, Hasinger, Priedhorsky 
& Middleditch 1989, Hasinger et al. 1990, Kuulkers et al. 
1994a, Kuulkers 1995, Kuulkers & van der Klis 1995, 1996, 
Kuulkers, van der Klis & Vaughan 1996) or magnetic field 
strength of the neutron star (Psaltis, Lamb & Miller 1995)0 
The Z sources CygX-2, GX5-1 and GX340-hO form one 
group, which we shall call the Cyg-like sources, while ScoX- 
1, GX 17-1-2 and GX 349-1-2 form another group, which we 
shall call the Sco-like sources. The Cyg-like sources have 
been proposed to have either a higher inclination, or a higher 
magnetic field, than the Sco-like sources. 

The only direct evidence that the compact stars in Z 
sources are neutron stars is the observation of X-ray bursts 
from GX 17+2 (Tawara et al. 1984, Kahn & Grindlay 1984, 
Sztajno et al. 1986) and CygX-2 (Kahn & Grindlay 1984, 
Kuulkers, van der Klis & Van Paradijs 1995, Wijnands et al. 
1997, Smale et al. 1996). The irregular occurrence of these 
X-ray bursts, and the lack of bursts in other Z sources has 
been attributed to the high mass accretion rate M in these 
systems (see the review by Lewin, Van Paradijs & Taam 
1993). 

Recently, the EXOSAT data of three Z sources have 
been analyzed in a detailed and homogeneous way (GX5-1: 
Kuulkers et al. 1994a, GX 340-1-0: Kuulkers & van der Klis 
1996, ScoX-1: Dieters & van der Klis 1996). The present 
work reports the results of a similar comprehensive study 
of the broad-band X-ray spectral and timing behaviour of 

t If it is due to orbital inclination the Z sources would be intrin- 
sically the same, but appear differently; if it is due to magnetic 
field, then they would be intrinsically different. 



GX 17-f2. Parts of the EXOSAT daXa. on this source have al- 
ready been reported by Langmeier et al. (1986), Langmeier, 
Hasinger & Triimper (1990), Sztajno et al. (1986), Stella, 
Parmar & White (1987), Schulz, Hasinger & Triimper (1989) 
and Hasinger & van der Klis (1989). In our study we dis- 
cuss the overall broad-band spectral and timing behaviour 
within the framework of the magnetospheric model for Z 
sources and of the idea that GX 17-1-2 is a Z source with a 
relatively low inclination. We also report previously unseen 
bursts, and place all EXOSAT hursts, within the Z-diagram 
(a preliminary report of these bursts was given by Kuulkers 
et al. 1994b). 



2 OBSERVATIONS 

We analysed all available EXOSAT medimn energy experi- 
ment (ME) argon data (Turner, Smith & Zimmerman 1981, 
White & Peacock 1988) of GX 17-1-2. In Table 1 we present 
the observation log of GX 17-1-2. The ME had eight detec- 
tors, each of which consisted of two layers, an argon filled 
chamber in front of a xenon filled chambei]^. Data were 
obtained in different observation modes, with good spec- 
tral and low time resolution (HER2, HER3, HER4, HER5; 
the different HER modes gave difi'erent options to transmit 
data from individual detectors or combined from more detec- 
tors), or with high time and no spectral resolution (HTR3, 
HTR4, HTR5). Alternatively, one could obtain data in two 
(e.g. 1985 day 258/259) or four energy bands (e.g. 1986 day 
093/094) at a high time resolution (HER7). A HER mode 
could be run simultaneously with one of the HTR modes. 
Either one half of the array of eight detectors ("half 1" [HI] 
or "half 2" [H2]) was pointed at the source while the other 
half monitored the background, or the whole array (WA) 
was pointed at the source (see Table 1). 

The HER modes were processed by the EXOSAT on- 
board computer (OBC) , which produced large dead-time ef- 
fects. In our analysis of CDs and HIDs (see below) we de- 
termine the dead-time factor from the qualified event rates 
(collected each 32s for each detector), which are (almost) un- 
affected by dead time. Of the HTR modes, only the HTR4 
mode is processed by the OBC. The HTR4 mode provides 
the count rate in only one (chosen) energy band at the high- 
est time resolution available with EXOSAT (~0.25 ms). The 
HTR4 mode observations on 1985 day 258/259 and on 1986 
day 093/094 provided the count rate in the 1.4-11.4 keV 
band and in the 1.2-19.9 keV band, respectively. 

On 1986 February 13th (Thursday, day 044) prob- 
lems due to atmospheric drag caused by orbital decay 
started with the EXOSAT aXtiiuAe and orbit control system 
(AOCS), see White (1986). Between March 27th (day 086) 
and April 9th (day 099) stable pointing was regained within 
5 arcmin. This resulted in a varying collimator response. In 
principle one could use the star tracker attitude control data 
(collected every 8 or 16 s) to correct for this. However, due to 
problems during the 1986 day 093/094 observations with the 
star trackers and gyros this was not possible for this obser- 
vation period. We therefore used these data only to produce 

t Early on 1985 day 232 one of the eight detectors ceased to 
operate. 
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a CD and not for the HIDs. After April 9th (day 099) stable 
pointing was lost again, which resulted in a complete loss of 
the satellite. On May 9th (day 129) EXOS AT re-entereA the 
Earth's atmosphere (White 1986). 



3 ANALYSIS 

For the spectral analysis we use CDs and HIDs. These di- 
agrams are created using either three or four energy bands 
(e.g. Hasinger & van der Klis 1989). The observations during 
1986 day 093/094 provided four energy bands, so we decided 
to analyse all the data using energy boundaries correspond- 
ing to these bands. We found that three energy bands (i.e. 
combining two of the four bands) gave the best separation of 
the branches of GX 174-2 in the CD: 1.2-4.7 keV, 4.7-6.6 keV 
and 6.6-19.9 keV. In the CD and HIDs the soft colour is de- 
fined as the ratio of the count rate in the second energy 
band to the count rate in the first energy band, while the 
hard colour is defined as the ratio of the count rate in the 
third energy band to that in the second energy band. The 
'intensity' used in the HIDs is the count rate per cm^ in the 
1.2-19.9 keV band, corrected for dead time and collimator 
response (as redetermined by Kuulkers et al. 1994a). The 
count rates were corrected for background. All points in the 
CD and HIDs are 200 s averages. 

In order to investigate the fast timing behaviour of 
GX 17-1-2 we performed fast Fourier transforms (FFTs) on 
successive 16 s and 512 s blocks of data of all the high-time 
resolution argon data. Power spectra of blocks of 512 s of 
data were used to study primarily the low-frequency be- 
haviour <lHz. In order to obtain a more accurate binning 
in position in the Z, we used the power spectra of 16 s blocks 
of data. The latter power spectra were used to study primar- 
ily the high frequency (>1 Hz) behaviour. The power spectra 
were grouped according to position in the Z and then aver- 
aged. 

When analysing power spectra of Z sources, one encoun- 
ters several power spectral components (Hasinger & van der 
Klis 1989): a constant level due to counting noise (the so- 
called Poisson level) modified by dead-time processes, three 
"noise" components, VLFN, LFN and HFN, and two kinds 
of QPO, HBO and NBO/FBO. In some cases the HBO 
show a second harmonic, which is also represented by a 
Lorentzian. For completeness we give the functional shapes 
of these components (see also Hasinger & van der Klis 1989, 
Kuulkers et al. 1994a): 

The VLFN is given as: PvLFN(i^) = AyI'~°''^ , where u is 
the frequency, ay the power-law index, and Ay the normal- 
ization constant. 

The LFN and HFN are given as: PNoisc(i') = 
A-mv~°'^ e~"^"^ , where q:n is the power-law index, the 
cut-off frequency, and the normalization constant. In the 
Z sources on was usually found to be consistent with zero 
for the HFN component (Hasinger & van der Klis 1989, see 
also Dieters & van der Klis 1996), and we fixed the parame- 
ter at this value. In Tables 2, 3a and 3b we use the subscript 
L and H for the LFN and HFN components, respectively. 

The HBO, their harmonic, the NBO and the FBO are 
broad peaks with centroid frequencies Vc- These peaks are 
described by Lorentzians: Pqpo(i^) ~ Aq\{v — vof' + 



(Ai^q/2)^] ^, where A/zq is the full width at half maximum 
(FWHM) of the QPO and Aq a normalization constant. 

The strengths of the various noise components are expressed 
in terms of the fractional rms amplitudes of the correspond- 
ing fluctuations in the time series. They are calculated by 
integrating their contributions to the power spectra over cer- 
tain frequency ranges, as determined from fits of the func- 
tional shapes given above to the power spectra. The VLFN 
was integrated in the 512 s power spectra over the range 
0.001 to IHz, and in the 16 s power spectra over the range 
0.01 to 1 Hz. For the LFN and HFN we integrated over the 
range 0.01 to 100 Hz. Note that in the 16 s power spectra 
we extrapolated the power spectral shapes from ~0.06 Hz 
down to 0.01 Hz. The HER7 and HTR4 data suffer from 
large variable dead time effects. We corrected the fractional 
rms amplitudes of the power spectral components for these 
effects (see van der Klis 1989, Kuulkers et al. 1994a). 

The Poisson levels of the power spectra from non-OBC 
processed data, the HTR3 and HTR5 modes, were deter- 
mined using the known relation between Poisson level and 
observed raw count rate for these modes (Van der Klis 1989, 
Berger & van der Klis 1994). In power spectra from OBC- 
processed data, the HER7 and HTR4 modes, the OBC dead 
time strongly affects the Poisson level (see e.g. Kuulkers et 
al. 1994a). Moreover, the power spectrum of the Poisson 
noise in these power spectral data can no longer be con- 
sidered constant at frequencies higher than ~120Hz; it in- 
creases to higher frequencies (Tennant 1987). Most HER7 
mode data had a Nyquist frequency of 128 Hz or lower, so 
this only had a significant effect in the HTR4 mode, which 
had a Nyquist frequency of 2048 Hz. We used the Poisson 
level as a free parameter in the fits to the HER7 power spec- 
tral data. In the case of HTR4 we first binned the data to a 
time resolution of ~4 ms, and used the Poisson level as a free 
parameter in the power spectral fits. We refer to Appendix 
A for a more detailed discussion of dead time effects. 

Whenever a power spectral component was not clearly 
present in the averaged power spectra, we determined a la 
upper limit to the rms amplitudes of this component by 
fixing its other parameters (while keeping the parameters of 
the other observed components free). The values of the fixed 
parameters were determined by evaluating the observed be- 
haviour (using only the HTR-mode data) of the correspond- 
ing components when they were clearly detected, and are 
discussed below. 

VLFN is present in all branches. Since there is no simple 
dependence of the VLFN power-law index as a function of 
the position in the Z (Section 4.2.2) we fixed it to its average 
observed value when determining an upper limit, i.e. 1.5. 

LFN is observed to be only present in the HB and the 
upper part of the NB. Upper limits where therefore only de- 
termined for Sz (see Section 4.2.1) up to ~1.35. No clear de- 
pendence of the LFN power-law index and cut-off frequency 
on position in the Z is found (Section 4.2.3); the fixed val- 
ues when upper limits were determined of these parameters 
were, therefore, —0.5 and 2 Hz, respectively. 

HFN is also present in the whole Z, but no clear depen- 
dence of the HFN cut-off frequency as a function of the posi- 
tion in the Z can be found (Section 4.2.4). In cases when an 
upper limit was determined the cut-off frequency was fixed 
at ~77. 
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The HBO centroid frequency clearly depends on position 
in the HB (Section 4.2.5, see also Penninx et al. 1990). The 
power spectral fit results from HTR3-mode data on 1985 
day 258/259 (Table 3a) indicate the following linear rela- 
tion between the centroid frequency and Sz'- i^hbo (Hz) = 
24.0 + 12.5 Sz- This relation was used when an upper limit 
on the HBO was determined; for the HBO FWHM no clear 
dependence on the position in the Z is found (Section 4.2.5) 
and we then fixed it at 3. 

A clear HBO harmonic was only detected on a few occa- 
sions (see Table 4, later). Whenever the fundamental HBO 
was clearly detected we determined an upper limit on its har- 
monic. In these cases the harmonic centroid frequency and 
harmonic FWHM were fixed to twice the value of the fun- 
damental frequency and fundamental FWHM, respectively, 
whenever the source was in the HB. Recently, QPO between 
50 and 62 Hz were found in the upper and middle part of the 
NB of GX 17-1-2 (Wijnands et al. 1996a). We regard this as 
the HBO harmonic (see Section 6.2.2). Between Sz values of 
1-1.75 we determined rms amplitude upper limits for such 
QPO using fixed values for the centroid frequency of 60 Hz 
and FWHM of 5 Hz. 

Upper limits to the NBO rms amplitude were determined 
in the entire NB. Average values of 3 Hz and 7 Hz were then 
used for the fixed centroid frequency and FWHM, respec- 
tively. 

Upper limits to the FBO rms amplitude were determined 
in the entire FB. We then fixed its centroid frequency and 
FWHM to the observed value, i.e. 16 Hz and 7 Hz, respec- 
tively. 



4 SPECTRAL AND TIMING BEHAVIOUR 

4.1 Broad-band spectral analysis 

The CD of all HER observations after 1983 is shown in 
Fig. la. In Figs, lb and c we display the corresponding HIDs, 
with the 1986 data, which had pointing problems, left out. 
The CD and HIDs of the 1983 day 215/216 data can be 
found in Figs. 2a-c. The data points in these figures have 
been connected to show more clearly the presence of the FB. 
It appears that the 1983 Z track is at higher soft and hard 
colours than that of the other observations. The overall in- 
tensities are comparable. CDs of other sources show similar 
shifts between data for a certain period in 1983 and later 
data. This was observed in the Z sources GX5-1, CygX- 
2, GX 349-1-2, and the atoll sources 4U 1636-53 and SerX-1 
(see Kuulkers et al. 1996, and references therein). Since all 
sources show the same effect, the shift is probably instru- 
mental. 

AU data on GX 17+2 obtained after 1983 faU within 
a single Z track in the CD. We can exclude shifts in the 
soft and hard colours at the level of ~2% (determined by 
examining the width of the branches) . We therefore conclude 
that the Z track of GX 17-1-2 in the CD is stable to within 
~2 per cent and does not show evidence for secular motion 
in the CD for observations separated by ~1 day up to ~1.5 
years. 

The HIDs of GX 17-1-2 seem to be also quite stable. 
However, one data set (1985 day 259) is placed at ~8 per 
cent higher intensities than the other data. Since instrumen- 
tal changes between the 1984 day 249/250 and the 1985 day 



258/259 data can only be expected to contribute at most up 
to 3 per cent in intensity (see Kuulkers & van der Klis 1996) 
we think this shift in intensity is probably intrinsic to the 
source. More observations are needed to verify if this effect 
is real and how it relates to other source properties. 

As can be seen in Figs, la-c and 2a-c, GX 17-1-2 was ob- 
served in all three spectral branches (see also Table 1). The 
HB is clearly not horizontal in the CD, but nearly vertical, as 
noted previously (e.g. Hasinger & van der Klis 1989); in the 
other diagrams it has different orientations. Due to an unfor- 
tunate choice of HER7 energy boundaries when the 1986 ob- 
servations were performed, the FB is positioned close to the 
NB (Fig. 4a). This also caused the FB to be oriented above 
the NB in Fig. 2b, while it falls along the NB in Fig. 2c. 



4.2 Power spectral analysis 
4-2.1 Introduction 

We studied the power spectral components as a function of 
position along the Z. We used a one-dimensional parameter, 
Sz to measure curve length along the Z (Hertz et al. 1992, 
Dieters & van der Klis 1996). The scale in the CD was fixed 
by assigning Sz values of 1 and 2 for the HB/NB apex and 
the NB/FB apex, respectively, in the CD of Fig. la. The 
same scaling was also used for the 1983 day 215/216 obser- 
vations. We note that during the first part of the 1985 day 
258/259 observation (16:15-03:13 UT) we could not estimate 
Sz from the CD, since data were only available in one broad 
low energy band. However, since the position in the HB is 
dependent on intensity we could determine Sz indirectly. 
For the second part of the 1985 day 258/259 observation we 
had sufficient data to determine the dependence of intensity 
on Sz in the relevant channel. We used this relation to de- 
termine Sz values in the first part of the observations. Using 
the above defined scaling we find that Sz ranges from ^0 in 
the upper part of the HB, to ~4.5 in the upper part of the 
FB. In Figs. 3a-c we show average power spectra from the 
different branches (HB, NB and FB). The various spectral 
components are indicated. 

In Tables 2, 3a and 3b we give the results of the power 
spectral fits as a function of Sz for the 512 s and 16 s data 
segments, respectively. In Table 4 we give the results for the 
HBO harmonic (see Section 4.2.5). Upper limits (la) have 
been determined on the rms amplitudes whenever a compo- 
nent was not clearly present in the average power spectra. 
The range in Sz where these upper limits where determined 
and the values to which the other parameters of that par- 
ticular component were fixed are described in Section 3. 

In Figs. 4a-l most of the power spectral parameters 
have been plotted vs. Sz. In the next sections we discuss 
the results for each power spectral component. The 1985 
day 258/259 HER7 and HTR4 data are not used in our 
discussion since they have different energy bands than the 
other data (see Table 1); only for the discussion of the HBO 
FWHM and frequency we include the HTR4 data, since 
these parameters are approximately independent of the en- 
ergy (Penninx et al. 1990, Penninx et al. 1991, Lewin et al. 
1992). 
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4-2.2 Very-low-frequency noise (VLFN) 

The VLFN was detected in all spectral branches, consistent 
with previous observations. The VLFN rms and its power 
law index are plotted vs. Sz in Figs. 4a and b, and Figs. 4c 
and d, for the 512 s FFTs and 16 s FFTs, respectively. We 
did not include the 1986 day 093/094 data, which had VLFN 
rms values between 5 and 6 per cent (Table 2), while in the 
other observations in the lower part of the NB the rms is 
around 1 per cent. We attribute this to the intensity varia- 
tions arising from the changing collimator response on time 
scales longer than ~100s, which causes leakage of power 
from lower frequencies into the passbands of the power den- 
sity estimators with nominal frequencies in the VLFN fre- 
quency range (see Deeter 1984). The power law index during 
the 1986 093/094 observations is about 2, consistent with 
low frequency leakage. 

The changing collimator response only affects frequen- 
cies below 0.01 Hz. This is clear from Figs. 4c and d, where 
the 1986 day 093/094 VLFN components, here measured 
only down to 0.01 Hz, fall along the other data points in 
both figures. 

The VLFN rms is lowest (~0.5 per cent) in the upper 
part of the HB and the upper part of the NB. It shows a 
local maximum in the middle part of the HB of ~1 per cent. 
From the upper part of the NB to the lower part of the 
NB the rms increases to ~1.5 per cent. From the 16 s FFTs 
we infer that the VLFN strength continues to rise from the 
lower part of the FB to the upper part of the FB, up to at 
least 3 per cent. 

In the lower part of the NB the power law index has 
values of ^--^1.5. In the upper part of the FB the power law 
index is about 2, as expected from leakage effects from slow 
(source) intensity variations (see above). 

4.2.3 Low-frequency noise (LFN) 

The LFN is peaked (see e.g. Fig. 3a) as is clear from the 
negative LFN power law index values between —1.5 and 0, 
with cut-off frequencies between 1 and 3 Hz. In Figs. 4e, f 
and g we show the LFN rms, power law index and cut-off 
frequency, respectively, as a function of Sz. The rms of the 
LFN decreases from ~3 per cent in the upper part of the 
HB to ~1.5 per cent in the upper part of the NB (see also 
Tables 2 and 3a). 

The low rms and the large uncertainty in the LFN cut- 
off frequency in the upper part of the NB indicate that the 
LFN vanishes in this part of the NB. In the rest of the NB 
and in the FB no evidence for LFN was found. 

4.2.4 High-frequency noise (HFN) 

The HFN is also present in all branches, see Figs. 4h and i 
for the 512 s FFTs, and Table 3a and 3b for the 16 s FFTs. 
Like the VLFN, the HFN is not correlated with any of the 
other power spectral components. The rms strength clearly 
decreases from the upper part of the HB to the middle part 
of the NB from ~5.3 per cent to ~2.5 per cent. From there 
onwards it remains between 2 and 3 per cent. 

The cut-off frequency is not well determined in many of 
the power spectra, mainly due to the relatively low Nyquist 
frequencies (64 Hz) in the 1984 day 249/250 and 250/251 



data. Whenever it could be determined we found values be- 
tween 10 and 50 Hz. 



4.2.5 QPO: HBO, NBO and FBO 

HBO were only found in the 1985 day 258 data in the upper 
part of the HB (Tables 2, 3a and 3b, Fig. 3a). There are 
indications of a very weak HBO component in the expected 
frequency range in the 1984 day 249/250 data, but this com- 
ponent is not significant. The QPO frequency increases from 
~23 to ~28 Hz from the upper part of the HB to the middle 
part of the HB, while the FWHM and rms decrease from ~7 
to ~2Hz and ~2.5 per cent to ~1 per cent, respectively. The 
HBO disappear below ~1 per cent rms in the lower part of 
the HB. 

A harmonic to the HBO peak with a frequency consis- 
tent with twice the HBO frequency is seen at the lowest Sz 
values (see Table 4, Fig. 3a). The harmonic has a FWHM 
between 10 and 20 Hz and an rms of about 2.5 per cent. The 
harmonic is about as strong as the HBO itself; the data do 
not allow to distinguish between the harmonic peak being 
the same width, or twice as wide, as the main peak (see also 
Table 4). 

NBO with frequencies between 5 and 7 Hz were only de- 
tected in the middle and lower part of the NB (Fig. 3b), with 
an rms between 1.5 and 2.5 per cent, and FWHM between 
2 and 6 Hz (Figs. 4j-l). 

In the 1986 day 093/094 data we found evidence for FB 
QPO (FBO) in the lower part of the FB (Fig. 3c), with a 
centroid frequency of ~16Hz, FWHM of ~7Hz and rms of 
~2 per cent. 



5 BURST ANALYSIS 
5.1 New bursts 

In Figs. 5a-e we show the light curves of all observations 
(except 1986 day 093/094) at 5 s time resolution. In Figs. 5c 
and d the smooth overall variations are interrupted by two 
bursts, which were already reported by Sztajno et al. (1986). 
In Figs. 5a-e we indicated which spectral branch the source 
was in at each instant. 

As noted in the previous section, the pointing during 
observation 1986 day 093/094 was only stable to within 
5 arcmin. In Fig. 6 we display this data set (Fig. 6a, 5 s 
time resolution) together with a soft-colour curve (Fig. 6b, 
200 s time resolution) . Although variations due to the vary- 
ing collimator response are visible, one can also see that the 
variations in the FB are more irregular than in the NB. 

Two sharp peaks seem to disturb the relatively smooth 
variations in the light curve of the 1986 day 093/094 data 
(~36 000 and ~104000s after the start of the observations). 
We investigated these peaks and found that they are X-ray 
bursts which had been previously overlooked. The house- 
keeping data ruled out a solar particle origin (see also e.g. 
Kuulkers et al. 1995). In Figs. 7a and b we display light 
curves of these bursts at a higher time resolution (0.5 s); 
plotted is the total HTR5 count rate, corrected for dead 
time and background. The second burst straddled two EX- 
OS AT ohserv&tion segments and is therefore interrupted for 



© 1996 RAS, MNRAS 000, 



6 Kuulkers et al. 



~200 s. The housekeeping monitoring during this time in- 
terval still provided 32 s resolution qualified event rate data 
with which we overlayed the HTR5 light curve of this burst. 

The bursts are not related to the changes in the pointing 
of the satellite. The events rise within 1 s, which is much 
faster than the pointing variations (typically ~100s), and 
show clear hardness variations, which can not be caused by 
pointing variations. In Figs. 7c and d we display the 4.7- 
19.9 keV/1.2-4.7 keV hardness ratio curves of the two bursts. 
The bursts display spectral softening during their decay. 

5.2 Burst properties 

The first of the two new bursts started at 1986 day 093 
14:04:57 (UT) and lasted for ~100 s, while the second started 
-19 hr later at 1986 day 094 9:05:40 (UT) and lasted ~300s. 
No other bursts occurred in between these two during times 
that the OBC was operational. We checked whether bursts 
occurred when the OBC was not operational using the 
housekeeping data, and found none. Since the collimator re- 
sponse changed on time scales of ~100 s during these bursts, 
it was not possible to accurately determine the standard X- 
ray burst parameters (see Lewin et al. 1993) a (ratio of the 
average persistent flux to the time-averaged flux emitted in 
the bursts) and 7 (ratio of the mean persistent pre-burst 
flux and the net peak burst flux). However, it is clear that a 
was very large (larger than 1000) and 7 > 1. Spectral data 
(HER2) were only available during the second burst, but be- 
cause the burst was interrupted and was observed with an 
uncertain, varying collimator response, we did not analyze 
these data. 

To investigate in which branch the source was when the 
1986 bursts and the bursts detected by Sztajno et al. (1986) 
occurred, we indicated the position of the soft and hard color 
values in the CD at the time the bursts occured with arrows 
in Fig. 1. We found that the smaU 1984 day 250 burst (I) 
occurred when GX 17-1-2 was in the upper part of the NB, 
while the other bursts (1985 day 232: II; 1986 day 093: III; 
1986 day 094: IV) occurred in the middle part of the NB or 
near the NB/FB apex. Although none of the bursts occurred 
in the HB, we note that GX 17-1-2 is only occasionally found 
in the HB (see Section 6.1). Moreover, using the total time 
spend in each branch as observed with EXOSAT, we found 
that the occurrences of the four EXOSAT burst are con- 
sistent with being independent of branch position. We also 
found no clear relation between burst duration and position 
in the Z. 

5.3 Pulsations during the bursts? 

We searched for periodicities in the bursts. We followed the 
same procedure for searching pulse periods as Vaughan et al. 
(1994). We transformed 16 s, 256 s and 64 s data segments for 
the 1984, 1985 and 1986 bursts, respectively. We also aver- 
aged the power spectra of 64 s data segments of the two 1986 
bursts, and of these bursts together with the 1985 burst. We 
found no evidence for a periodic signal above the 99 per cent 
confidence detection level. The 99 per cent confidence upper 
limits to the amplitude of a sinusoidal modulation are given 
in Table 5. The best limits are ~1 per cent for the 0-256 Hz 
frequency range, and ~1.5 per cent for the 256-512 Hz fre- 
quency range. 



6 DISCUSSION 

6.1 Two groups of Z sources 

GX 17+2 and the other Z sources form a group with similar 
behaviour, i.e. three branches with correlated fast timing 
behaviour (Hasinger & van der Klis 1989). However, it has 
been shown that GX 17-1-2 shares properties with ScoX-1 
and GX 349+2, which are distinct from GX5-1, CygX-2 
and GX 340+0 (Hasinger & van der Klis 1989, Penninx et al. 
1991, Kuulkers et al. 1994a). The latter sources, for example, 
show long-term variations in the position of the Z pattern 
in the CD by typically —7 per cent (which are recurrent in 
position and shape in the case of CygX-2, see Kuulkers et 
al. 1996) and dipping behaviour in the FB, while the former 
sources show a more stable Z pattern and flares in intensity 
in the FB. In this paper, we showed that in the CD the 
position of the Z track of GX 17+2 is stable to within 2 per 
cent, which is in accordance with this. 

It has been suggested that these differences are caused 
by a difference in orbital inclination angle (Kuulkers et 
al. 1996, and references therein). Optical observations of 
CygX-2 (Cowley, Crampton & Hutchings 1979) and ScoX-1 
(Crampton et al. 1976) indicate orbital inclinations of 65- 
75° and 15-40°, respectively. The fact that intensity de- 
creases in the FB of the Cyg-like sources is interpreted as 
obscuration of the central emission by the inner disk when it 
swells up as the luminosity approaches the Eddington limit. 
Orbital inclination also plays a key role in suggested ex- 
planations for: (i) the secular variations in the shape and 
position in the CD and HID of the Z pattern of the Cyg-like 
sources (e.g. Kuulkers et al. 1994a, 1996), (ii) the shortness 
of the FB in GX 5-1 (Kuulkers & van der Klis 1995) and (iii) 
the pecuhar QPO found in the FB of CygX-2 (Kuulkers & 
van der Klis 1995). For example, long term variations may 
be caused by changes in our view of the inner regions due to 
changes in the geometry of structures near the orbital plane 
(such as accretion stream and outer disk, see Kuulkers et al. 
1994a, 1996, see also Wijnands, Kuulkers & Smale 1996b). 

Several other properties (described below) in Z sources 
also appear to correlate with the division into two groups, 
but have not (yet) been interpreted in terms of orbital incli- 
nation differences. 

Peaked LFN has been observed in GX17+2 and ScoX- 

1 (see e.g. Hasinger & van der Klis 1989), while all three 
Cyg-like sources show non-peaked LFN. Atoll sources also 
show peaked and non-peaked noise (the so-called atoll source 
HEN). It has been suggested that atoll source HFN is similar 
to Z source LFN (van der Klis 1994). The different HFN 
shapes in atoll sources could perhaps be a manifestation of 
similar effects to those seen in Z source LFN. 

A vertical (short) HB in the Sco-like sources has been 
seen, versus a horizontal HB (with an upward bend at the 
lowest M) in the Cyg-like sources. 

The following property also appears different between the 
two groups, but seems hard to fit in with the orbital inclina- 
tion framework: GX 17+2 and ScoX-1 are only occasionally 
found in the HB, and GX 349+2 never, while GX 5-1, Cyg X- 

2 and GX 340+0 are frequently found in the HB, but only 
occasionally in the FB. Moreover, the Sco-like sources spend 
a much longer time scale in the FB (>12hr) than the Cyg- 
like sources (<2hr). Since there is a good evidence that M 
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increases from the NB to the FB this imphes that, on aver- 
age, the Sco-Hke sources have a higher M than the Cyg-hke 
sources. It is hard to understand how this could have its ori- 
gin in the inchnation at which we view the sources, unless 
viewing geometry affects the M level at which branch transi- 
tions occur in the CD and HIDs. Careful study of fast timing 
properties as a function of Sz can in fact show whether this 
is the case, but the available data are not sufficient to make 
a strong case either way. 

Recently, Psaltis et al. (1995) suggested the magnetic 
field strength of the neutron star to be the origin of the 
differences in the shape of the Z track in the Cyg-like and 
Sco-like sources (the magnetic field being lower in the lat- 
ter sources). From their radiative transfer calculations they 
infer that subtle differences in the magnetic fields are able 
to explain the observed differences in the shapes of the Z 
pattern. It is not clear whether a difference in neutron star 
magnetic fields can explain the shifts in the Z track position, 
the apparent difference in mean M between the two groups, 
or the differences in LFN shape, but it does seem likely that 
differences in HBO properties between the two groups can 
be explained by it. 



6.2 Timing behaviour 

6.2.1 Very-low-frequency noise 

The VLFN was detected in aU branches. Near the HB/NB 
vertex it was weakest (~0.5 per cent rms), while in the FB 
it was strongest (rms of up to at least ~3 per cent rms). We 
found power law indices between ~0.8 and ~2. 

The VLFN of Z sources has been investigated in detail 
in several other studies as a function of Sz (GX5-1: Ku- 
ulkers et al. 1994a, EXOS AT argon; ScoX-1: Dieters & van 
der Klis 1996, EXOSAT xenon; Hertz et al. 1992, Ginga; 
GX 340+0: Kuulkers & van der Klis 1996) or as a function 
of rank number (GX5-1: Lewin et al. 1992, Gmga; CygX-2; 
Hasinger et al. 1990, Wijnands et al. 1996c, Gmga; Hasinger 
1993 [private communication], Ginga and EXOSAT argon; 
GX 340-1-0: Van Paradijs et al. 1988a, EXOSAT &vgon; Pen- 
ninx et al. 1991, Gmga; GX 17-1-2: Penninx et al. 1990, 
Ginga). Hasinger & van der Klis (1989) investigated EX- 
OSAT data of the six Z sources as a function of HB, NB 
and FB. Comparing the EXOSAT argon and the Gmga re- 
sults (both instruments provided roughly the same energy 
range), we find that the VLFN is lowest in the upper part of 
the NB, near the HB/NB vertex. The minimum strength is 
similar for all the sources, i.e. near 0.5 per cent rms. The in- 
crease in VLFN from the upper part of the NB to the lower 
part of the FB is strongest in GX5-1 (Kuulkers et al. 1994a) 
and CygX-2 (Hasinger 1993, private communication). The 
VLFN rms near the NB/FB vertex is ^4 per cent in GX5- 
1 (Kuulkers et al. 1994a), while it is about 2 per cent for 
GX 340-1-0 (Van Paradijs et al. 1988a, Penninx et al. 1991), 
GXn+2 (this paper) and ScoX-1 (Hertz et al. 1992). The 
sources GX 5-1 and GX 17-1-2 show evidence for a maximum 
of ~1 per cent in the middle part of the HB (Kuulkers et al. 
1994a and this paper, respectively). As shown by Dieters & 
van der Klis (1996), the VLFN is stronger in the 5-35 keV 
xenon data (see also Stella et al. 1987), i.e. at higher ener- 
gies. 

The power law indices for GX 340-1-0 (Penninx et al. 



1991), GX 17-1-2 (Penninx et al. 1990, this paper) and ScoX- 
1 (Hertz et al. 1992, Dieters & van der Klis 1996) show more 
or less the same behaviour as a function of Sz. The index 
is low (;$1) in the HB and upper part of the NB, while it 
increases to ~2 near the NB/FB vertex. In the FB VLFN 
power law indices around 2 have been found (note that the 
true values could be even higher, due to the leakage of power 
in the VLFN frequency range). GX 349-1-2 (Hasinger & van 
der Klis 1989) shows behaviour consistent with this. GX 5-1 
(Kuulkers et al. 1994a) shows a clearly opposite behaviour. 
The power law index remains at about 2 along the HB to 
the middle part of the NB, and then drops to ~1.4 in the 
lower part of the FB. 

In the model of Lamb (1989) the increase of VLFN from 
the lower part of the NB into the FB is explained by instabil- 
ities in the accretion flow. When M reaches the Eddington 
limit radiation pressure becomes important. Accretion onto 
the neutron star becomes highly unstable at this stage, and 
therefore increases the VLFN strength. 

6.2.2 Horizontal branch QPO 

In the upper/middle part of the HB we found HBO with 
frequencies between 24 and 28 Hz, FWHM between 1 and 
8 Hz, and fractional rms amplitudes between 1.5 and 2.5 per 
cent. The HBO decreased rapidly in strength from the up- 
per part of the HB to the middle part of the NB. We found 
evidence for a second harmonic, with a frequency consistent 
with twice the HBO frequency, and with similar strength 
as the first harmonic. These results are similar to those re- 
ported by Stella et al. (1987) and Langmeier et al. (1990) 
for the EXOSAT data. The range of the HBO frequencies 
falls within the range (18-30 Hz from upper to lower part of 
the HB) found by Penninx et al. (1990) in Ginga data. 

More recently, variable- frequency QPO have been found 
on the NB of GX17-f2 (Wijnands et al. 1996a). These au- 
thors interpreted these QPO as being (fundamental) HBO 
which persists on the NB. However, if we extrapolate the 
(linear) relation between HBO fundamental centroid fre- 
quency and S'z, using our HTR-mode EXOSAT data, we 
find a frequency of ~36.5Hz at the HB/NB vertex (Sz^l), 
which is near the maximum value found in the lower part of 
the HB in the Ginga data (see above). A similar conclusion 
can be derived from the Ginga power spectral parameters for 
the HBO (Penninx et al. 1990). Extrapolating the (linear) 
relation between HBO fundamental centroid frequency and 
either Ginga MPC3 colour or Ginga PC colour to ~60 Hz 
give HB/NB vertex colour values which are way off from 
the expected place of the vertex in the Gmga CD and HID. 
We therefore interpret the variable-frequency QPO as be- 
ing a harmonic of the HBO, unless the HBO fundamental 
centroid frequency changes abruptly (by 15-20 Hz) near the 
HB/NB vertex. The linear change of the HBO harmonic 
centroid frequency (see Penninx et al. 1990) as a function 
of either Ginga MPC3 colour or Gmga PC colour, extrapo- 
lated to the expected HB/NB vertex in the Ginga CD and 
HID are consistent with this interpretation. 

The other Z sources in which HBO have been found, 
GX5-1 (Lewin et al. 1992, Kuulkers et al. 1994a), CygX-2 
(Hasinger 1987a, Hasinger et al. 1990, Wijnands et al. 1997) 
and GX 340-1-0 (Penninx et al. 1991, Kuulkers & van der 
Klis 1996) show HBO properties that are somewhat different 
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from those of GX 17+2. The maximum HBO frequencies are 
much higher in the former sources (up to 55 Hz, with rms 
amphtudes between 2 and 6 per cent) . In GX 5-1 a second 
harmonic with an rms amphtude roughly half that of the 
first harmonic was found (Lewin et al. 1992, Kuulkers et 
al. 1994a); similar harmonics are reported in the HBO of 
CygX-2 (Hasinger 1987a, Focke 1996, Wijnands et al. 1997) 
and GX 340+0 (Kuulkers & van der Klis 1996). 

The HBO are thought to be caused by magnetospheric 
gating (Alpar & Shaham 1985, Lamb et al. 1985, Ghosh 
& Lamb 1992), at the beat-frequency between the neutron 
star spin frequency (i^s) and the Kepler frequency at the 
inner disk edge (i^Ko)- This beat-frequency model requires 
inhomogeneities in the accretion flow, whose presence has as 
a natural consequence the production of LFN. The radius 
of the inner edge of the disk depends on M, the magnetic 
moment fi and the mass M of the compact object (Ghosh 
& Lamb 1992). For n-fold symmetry in the magnetic field 
pattern, this leads to (see Lamb et al. 1985, Ghosh & Lamb 
1992): 

^HBO = niCM'^ii^M'' - //Ns), (1) 

where n is the harmonic number, C is a constant, and a, 
P and K are exponents which depend on the particular in- 
ner disk model considered; a > 0, /3 < 0, whereas k can be 
either positive or negative. There are various ways to pro- 
duce the ^factor 2 lower HBO frequency range in GX 17+2 
as compared to the other Z sources: lower the overall mass- 
accretion rate, increase the magnetic field strength or the 
mass of the neutron star, either increase or decrease the spin 
period of the neutron star depending on the value of k, or 
lower the harmonic number n (if n>l). Of these possibilities, 
differences in M and M are unlikely. M is probably the same 
at the apexes for all Z sources (e.g. Lamb 1989, Hasinger et 
al. 1990). (But perhaps the inclination influences the posi- 
tion of the apexes in an unknown way, see also Section 5.2.) 
Unless K is large, neutron star masses are not likely to be 
sufficiently different to explain a factor of 2 in the maxi- 
mum HBO frequency. A closer match between i^ns and vko 
in GX 17+2 than in the other Z sources would of course also 
lead to a lower HBO frequency. However, in that situation 
the range in HBO frequency, which is determined by the 
range in vko, would still be expected to be similar to that in 
the other sources, whereas it is lower: ~12Hz in GX17+2 
versus ~35Hz in GX5-1 and CygX-2. 

Since the maximum HBO frequency in GX 17+2 is 
about a factor of two lower than the maximum HBO fre- 
quency in CygX-2, GX5-1 and GX 340+0, we now explore 
the possibility that the HBO frequency seen in GX 17+2 is 
dominated by the fundamental (first harmonic) of the field 
symmetry {n=l) and the HBO frequency seen in the lat- 
ter sources by the second harmonic (n=2). This could be 
true if in GX17+2 the radiation from one magnetic pole 
at the neutron star surface dominates that from the other 
magnetic pole, whereas in the other sources both poles con- 
tribute approximately the same amount of radiation. The 
less luminous magnetic pole in GX 17+2 is expected, if its 
luminosity is not negligible, to give rise to a second har- 
monic with centroid frequencies in a similar range as the 
HBO in the other sources, as actually observed in GX 17+2. 
The second harmonic seen in GX 5-1 and GX 340+0 is weak 
compared to that in GX 17+2 and in our interpretation due 



to deviations of the HBO signal from being sinusoidal. If the 
above described model is right, we predict the presence of 
a subharmonic (due to the difference between the emission 
from the two poles) in the power spectra of GX5-1, CygX- 
2 and GX 340+0, while it should be absent in the power 
spectra of GX 17+2. 

Differences in the radiation from the two magnetic poles 
could be caused by asymmetries in the field geometry, which 
lead to different rates of accretion onto, or different emis- 
sion characteristics from, the two poles. These asymmetries 
would also lead to a difference in the surface area of the two 
polar caps. Different polar cap areas have also been deduced 
from pulse analyses in several pulsars (e.g. Leahy 1991, Bulik 
et al. 1992). 

We note that in a situation with unequal polar cap ar- 
eas type I bursts might have a higher chance to occur, as 
on the larger of the two polar caps accretion per unit area 
would be less, and therefore on that pole the conditions for 
a type I burst would be more favourable than for two equal- 
area polar caps. The presence of bursts and the lower maxi- 
mum HBO frequency might therefore be explained together 
in this model. 

We shall for illustrative purposes consider one simple 
asymmetry in field geometry, namely that which arises when 
the (dipole) magnetic field center does not coincide with the 
neutron star center, but is shifted. In Figs. 8a and b we 
plotted a magnetic dipole field in which the neutron star 
center coincides with the magnetic field origin (Fig. 8a) and 
one in which the neutron star origin is somewhat shifted 
(Fig. 8b). In these figures we have indicated the polar caps 
as expected from the last closed field line in a simple dipole 
plus spherical magnetosphere geometry. It is known that de- 
viations from this geometry will occur, especially near the 
inner disk edge, because of the magnetic field of the accretion 
disk (e.g. Ghosh & Lamb 1979, Spruit & Taam 1990), but 
we ignore this here. We computed the areas of the two po- 
lar caps by assuming magnetospheric radii between 15 and 
100 km (which is roughly the range expected in Z sources, 
see Ghosh & Lamb [1992]), and offsets of the magnetosphere 
origin with respect to the neutron star origin between and 
5 km. The neutron star was assumed to have a radius of 
10 km. In Table 6 we give the results. As can be seen in this 
table, a slight offset of the magnetic dipole has a large effect 
on the polar cap areas. The ratio of the two polar cap ar- 
eas already exceeds 2 for offsets of only 1.5 km. For small 
magnetospheric radii, a small offset will cause the larger of 
the two polar caps to cover a considerable fraction of its 
hemisphere. 

When M increases along the Z, the magnetospheric ra- 
dius becomes smaller. This leads to a rapid increase in the 
polar cap areas, which might explain the fact that the four 
bursts observed with EXO SAT occurred in the NB and FB, 
and not in the (lower M) HB. However, somewhat depen- 
dent on the inner disk model (see Ghosh & Lamb 1992), 
accretion per unit area is usually still expected to increase 
when total M increases in our simple model, even on the 
larger of the two polar caps. 

The presence of multipole components in the magnetic 
field is another possibility to cause the required asymmetries 
between the accretion regions. As multipole effects increase 
very strongly with decreasing radius, these effects would 
be most prominent for small magnetospheric radius, i.e., at 
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higher M or for lower overall magnetic field strengths. An 
overall lower magnetic field in GX 17+2 (as suggested by 
Psaltis et al. 1995, see the discussion in Section 6.1) could 
also be the explanation of the occurrence of bursts in this 
source, and the increase of M from HB towards NB might 
explain why no bursts have been found in the HB. 

6.2.3 Flaring branch QPO 

On one occasion we found evidence for FBO with a fre- 
quency, FWIfM and rms amplitude of 16 Hz, 7 Hz and 2 
per cent, respectively. Similar FBO in GX 17+2 were re- 
ported earlier in Gmga data by Penninx et al. (1990) with 
frequencies up to 20 Hz and rms up to 3 per cent, but not 
from EXOSAT data. ScoX-1 is the only other Z source for 
which similar FBO have been reported (see Dieters & van 
der Klis 1996, and references therein). This kind of FBO 
merges smoothly with the NBO, and therefore probably has 
a common origin. FBO, and the NBO, are believed to be due 
to oscillations in the optical depth of the radial flow in the 
inner disk region (Lamb 1989, Fortner, Lamb & Miller 1989) 
or due to sound waves in the thick accretion disk (Hasinger 
1987b, Alpar et al. 1992). 

It was recently suggested that the strength of the NBO 
oscillations depends on the inclination at which we view the 
Z sources (Kuulkers & van der Klis 1995), stronger NBO 
occurring in the lower-inclination sources. The strength of 
the NBO in the lower NB and the fact that we observe FBO 
in GX 17+2 are in accordance with this. It is not clear how 
this difference could be explained in the Psaltis et al. (1995) 
model for the two groups of Z sources. 

6.2.4 High-frequency noise 

We found that in GX 17+2 the strength of the HFN de- 
creases as a function of Sz from the HB to the middle part 
of the NB from ~5.3 per cent to ~2.5 per cent. From the 
middle part of the NB to the FB the HFN is more or less con- 
stant at a level of ~2.5 per cent. This is different from what 
is seen in ScoX-1 (Hertz et al. 1992, Dieters & van der Klis 
1996), and probably CygX-2 (Hasinger et al. 1990, Hasinger 
1993, private communication). There the HFN is reported 
to decrease from the HB, to the NB, into the FB. None of 
these observations were obtained with the EXOSAT ME ar- 
gon detectors, as our observations. The results for ScoX-1 
and CygX-2 from Gmga (Hertz et al. 1992, Hasinger et al. 
1990, Hasinger 1993, private communication) indicated val- 
ues which were typically lower than ~1.5 per cent in the 
FB. It has recently been suggested that some of the HFN 
seen with the EXOSAT ME is instrumental, up to an rms 
strength of ~2.3 per cent (0.01-100 Hz) in the argon data 
(Berger & van der Klis 1994, 1997). This instrumental effect 
has an amplitude similar to the values found in the mid- 
dle part of the NB and FB of GX 17+2. We conclude that 
our data are therefore consistent with a decrease in HFN 
strength from the HB to the FB. 

6.3 Bursts in GX 17+2 

We have detected two bursts in GX 17+2 in the 1986 data. 
One burst was about 100 s long, the other (which occurred 



~19hr later) lasted ~300s. Previously, Sztajno et al. (1986) 
had found two bursts in the 1984 and 1985 data of EX- 
OSAT, which lasted ~10 and ~300s, respectively. Bursts 
in this source have also been observed with other X-ray in- 
struments: one with Einstein (~10s, Kahn & Grindlay 1984) 
and four with Hakucho (two which lasted for ~100 s and two 
which lasted for ~300s, Tawara et al. 1984). All bursts have 
large a values (q:>1000) and 7 larger than 1 (i.e., the net 
peak burst flux is smaller than the average persistent flux). 

We found that the EXOSAT bursts occurred in the NB 
and lower part of the FB, i.e. at near-Eddington mass accre- 
tion rates. None were found in the HB, i.e. when M is lowest, 
although we can not rule out the possibility that this is a 
coincidence. The only other Z source that has shown X-ray 
bursts is CygX-2 (Kahn & Grindlay 1984, Kuulkers et al. 
1995, Wijnands et al. 1997, Smale et al. 1996). The occur- 
rence of these bursts was found to be independent of branch 
position (Kuulkers et al. 1995, Wijnands et al. 1997). 

Bildsten (1993, 1995) suggested that rotating patches 
of nuclear-burning material on the neutron star surface dur- 
ing X-ray bursts may produce regular pulsations (see also 
Schoelkopf & Kelley 1991). We searched for such pulsations 
in all four EXOSAT huTsts, but found none, with upper lim- 
its to the modulation strength of ~1 per cent (if the signal 
is sinusoidal). A similar lack of pulsations was found for the 
persistent emission (Vaughan et al. 1994b) and in a large 
sample of X-ray bursts from other sources (Jongert & van 
der Klis 1996). 

From the fact that the two new EXOSAT bursts show 
spectral softening during their evolution, we conclude that, 
like the bursts reported by Sztajno et al. (1986), they are 
type I X-ray bursts caused by thermonuclear flashes on the 
neutron star surface (Hoffman, Marshall & Lewin 1978, see 
also Lewin et al. 1993 and references therein). The bursts 
in CygX-2, although much weaker than those in GX 17+2, 
are also thought to be due to thermonuclear flashes, since 
the burst properties resemble those of type I bursts (see 
Kuulkers et al. 1995, Wijnands et al. 1997). However, no 
evidence for spectral softening was found. A burst recently 
detected with RXTE in CygX-2 (Smale et al. 1996) is the 
first case of a bona fide type I burst in this source. The fact 
that two Z sources show type I bursts indicates that the 
compact star in systems of this type is a neutron star. 

High luminosity burst sources tend to show irregular 
bursting behaviour (Van Paradijs et al. 1988b). They show 
bursts after long and variable waiting times. This is most 
probably due to steady nuclear burning, which occurs at 
high mass accretion rate M. It is unlikely that M is lower in 
GX 17+2 and CygX-2 as compared with the other Z sources, 
as in the NB M is believed to be approximately Eddington 
in all Z sources (Lamb 1989). If bursts do apparently occur 
at these high levels of M (see e.g. Taam, Woosley & Lamb 
1996), then it is not clear why bursts are not seen in the other 
Z sources as well. From the presence of HBO in Z sources 
one derives that they contain a neutron star with a non- 
negligible magnetic field. In Section 6.2.2 we have discussed 
the possibility that the magnetic field strength might affect 
the probability to produce bursts in a Z source. 
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7 SUMMARY 

We found that the position of the Z pattern in the CD is 
not moving by more than a few percent between observa- 
tions separated by from 1 day up to ~1.5 years. This is 
in accordance with the division of the Z sources in two 
groups, in which the Cyg-like sources (GX5-1, CygX-2 and 
GX 340+0) show secular variations in their Z pattern in 
the CD, while the Sco-like sources (GX17+2, ScoX-1 and 
GX 349+2) do not. 

We have investigated the various power spectral compo- 
nents as a function of position in the Z, Sz, and placed their 
properties in the framework of the division of Z sources into 
two groups and discussed the hypothesis that this division 
is due to differences in viewing geometry. Many, but not all, 
source properties fit in this picture. 

Two new X-ray bursts were found in GX 17+2 lasting 
~100 and ~300 s, which came at an interval of ~19 hr. These 
bursts, and also those seen by Sztajno et al. (1986), occurred 
in the NB or in the lower part of the FB, i.e. presumably 
near Eddington accretion rates. No evidence for regular pul- 
sations was found in the bursts. 

We suggest that an asymmetric magnetic field, giving 
rise to different surface areas and emission characteristics of 
the magnetic pole areas at the neutron star surface may si- 
multaneously explain the occurrence of bursts in GX 17+2, 
the fact that they were not seen at the lowest accretion rates, 
the fact that the maximum HBO frequency in GX 17+2 is 
lower by about a factor two than the maximum HBO fre- 
quency in CygX-2, GX5-1 and GX 340+0 near the HB/NB 
apex, and the fact that the observed range in HBO frequency 
in GX 17+2 is about a factor of two smaller than that of 
CygX-2 and GX5-1. The asymmetry could for example be 
due to an overall lower field strength, allowing the multipole 
component of the magnetic field to play a larger role due to 
the smaller radius of the atmosphere. 
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APPENDIX A: DEAD TIME AND POISSON 
NOISE 

Dead time is a severe problem in X-ray astronomy (e.g. 
Lewin, Van Paradijs & van der Klis 1988, van der Klis 1989). 
Once an X-ray photon is detected, the detector remains 
"dead" for a certain or variable time. During this "dead 
time" the detector is unable to register new incoming X-ray 
photons. The various dead time processes onboard EXOSAT 
have been described by Andrews (1984), Andrews & Stella 
(1985) and Tennant (1987). When the detector is dead for a 
fixed amount of time after a photon has been detected, one 
speaks of a fixed dead time. In the case the detector can only 
register one photon within a given time (or "sample cycle"), 
the dead time is variable (i.e. depending on where the pho- 
ton is detected within the sample cycle), and one speaks in 
terms of variable dead time. Within the ME instrument (i.e. 
the detector, electronics and the OBC) both constant dead 
time (in non-OBC processed data, such as HTR3 and HTR5 
observation modes) and variable dead time (in HTR4 and 
all HER observation OBC modes) processes occur (Andrews 
1984, Andrews & Stella 1985, Tennant 1987, Berger & van 
der Klis 1994). 

Dead time processes have recognizable effects on power 
spectra (see Lewin et al. 1988, van der Klis 1989). Pure 
counting noise, in the absence of dead time, has a white 
power spectrum with, in the case of the Leahy et al. (1983) 
normalization, an average power of two. This white noise 
component is called Poisson noise and its level the Poisson 
level. The Poisson level is modified by the dead time process. 
The larger the dead time, the lower the Poisson level. Also 
the dead-time modified Poisson noise is not necessarily white 
any more. The effect of the fixed dead time of HTR3 and 



HTR5 on the power spectrum has been studied by Berger 
& van der Klis (1994, 1997). 

In the case of variable dead time, the Poisson level is 
approximated as follows (e.g. van der Klis 1989): 

Ppoiss = 2(1 - AtTsamplc), (Al) 

where Tsampic is the duration of the sample cycle, and the 
observed count rate. For the EXOSAT ME value of Tsampie, 
1/4096 s, this expression gives wrong results. For a so-called 
effective sampling time of 1/3569 s (Andrews & Stella 1985) 
an approximate match with the data is found. 

Recently we found (Kuulkers et al. 1994a) that for 
high observed count rates (1500-2300 ctss^^) in the OBC 
processed HER7 data the fitted Poisson levels deviated 
from that expected from Eq. Al. Since the HER7 data of 
GX 17+2 provided lower count rates (0-1500 cts s~^), we de- 
cided to analyse all four-channel HER7 data for the bright- 
est LMXBs, i.e. ScoX-1, GX5-1, CygX-2, GXn+2 and 
GX 3-1-1. The former four sources are all Z sources, while 
the latter is an atoll source (see e.g. Hasinger & van der 
Klis 1989). In Table Al, we give an observation log of the 
four channel HER7 argon observations of these sources. We 
selected the power spectra according to observed intensity 
and determined the Poisson level by fitting the power spec- 
tra with the various source components as described in Sec- 
tion 3, with the Poisson level as a free parameter. This was 
done for each source individually. The results are given in 
Fig. Al. In this figure we also plotted the Poisson levels as 
predicted from the linear relation Eq. Al. It is clear that 
the fitted Poisson levels deviate from those expected (espe- 
cially at high count rates), as was already found by Kuulkers 
et al. (1994a). It can be seen that the Poisson level devia- 
tions from the predicted linear relation in the various sources 
where they overlap are the same. The deviation is therefore 
not an effect related to variability in the sources themselves. 

As mentioned in Section 3, we rebinned the OBC pro- 
cessed HTR4 data to a time resolution of ~4ms, i.e. similar 
to most of the HER7 data. In the power spectral fits of the 
HTR4 data the Poisson level was also taken to be a free 
parameter. The results of these Poisson level fits are also 
plotted in Fig. Al. As can be seen in this figure, the HTR4 
data fall along the HER7 data points, as expected from the 
fact that these data are processed in a similar way, by the 
OBC, to the HER7 data. 

Although a Unear deviation was sufficient to describe 
the GX 5-1 data of Kuulkers et al. (1994a), they noted that 
the Poisson level could be better described as a quadratic 
function of observed intensity. Here we find that a linear 
fit is unacceptable, and that our data are described by a 
quadratic function. When fitting the data to such a function, 
the Poisson level as a function of observed intensity, n, is as 
follows: 

Pvoiss HER7 = 1.9959(19)-5.89(4)-10-V+3.04(13)-10-V',(A2) 

where we give the uncertainty in the last digit (s) in brackets. 
The fitted Poisson levels in the observed count rate range 
100-900 cts s~^ are below the expected Poisson levels. This is 
consistent with a previous study by Rutledge (1993 private 
communication) from HER7 data of the Rapid Burster. 

We note that it was possible to fit the following function 
to the data: 
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ftoiss HER7 = 2 ^1 + ^(-l)'i(MTsamplo)' j . (A3) 

We found that y^^d was ~1.0 for n=4. In that case Tsampio 
was found to be 1/3304 s. 



Figure Captions 

Figure 1: 

(a) EXOSAT X-ray colour-colour diagram and (b, c) 
hardness-intensity diagrams of GX 17-1-2. Soft colour is de- 
fined as the ratio of the count rates in the 4.7-6.6 keV and 
the 1.2-4.7 keV bands, hard colour as the ratio of the counts 
in the 6.6-19.9 keV and the 4.7-6.6 keV bands. The inten- 
sity (corrected for background, dead-time and collimator 
response) is defined as the count rate in the 1.2-19.9 keV 
band. Different symbols refer to different observation pe- 
riods, which are given in the lower right part of (a). Each 
point represents a 200 s average. Typical error bars are given 
in each frame. HB (horizontal branch) , NB (normal branch) 
and FB (fiaring branch) are indicated. In a we give the po- 
sition in the CD at the time the different bursts (1984 day 
250: I, 1985 day 232: II, 1986 day 093: III, 1986 day 094: IV) 
occurred (see Section 5.2). 

Figure 2: 

(a) Colour-colour diagram and (b, c) hardness-intensity di- 
agrams of GX17+2 of the 1983 data. The colours and in- 
tensity are as defined in Fig. 1. The points are connected 
to show more clearly the different branches, which are indi- 
cated. Each point represents a 200 s average. 

Figure 3: 

Three examples of average power spectra corresponding to 
different segments in Sz. The corresponding fits are shown 
for clarity, and the power spectral components are indi- 
cated, (a) Power spectra for data of 1985 day 258/259 when 
GX 17-1-2 was in the upper part of the HB near S'z~0.1. 
Both the HBO and its harmonic can be seen, as well as the 
peaked noise component LFN. (b) Power spectra for data 
of 1985 day 232/233 during the NB near Sz~1.8. The NBO 
can be seen near 7.6 Hz. (c) Power spectra for data in the 
lower part of the FB during 1986 day 093/094 near S'z~2.4. 
The FBO can be seen near 16 Hz. 

Figure 4: 

(a)-(l) Properties of the various power spectral components 
as a function of position on the "Z" , Sz . The different sym- 
bols refer to different observation periods as indicated in 
Fig. la. The 1983 day 215 data are indicated with Q. The 
boundaries between the spectral branches are indicated, (a) 
VLFN rms 512 s FFTs, (b) VLFN power law index 512 s 
FFTs, (c) VLFN rms 16 s FFTs, (d) VLFN power law in- 
dex 16 s FFTs, (e) LFN rms 16 s FFTs, (f) LFN power law 
index 16 s FFTs, (g) LFN cut-off frequency 16 s FFTs, (h) 
HFN rms 512 s FFTs, (i) HFN cut-off frequency 512 s FFTs, 
(j) NBO frequency, (k) NBO rms (1) NBO FWHM. 



Figure 5: 

EXOSAT ME light curves of 4 observation periods of 
GX17-f2 at a time resolution of 5 s. Intensity is the dead- 
time-, background- and collimator-response-corrected count 
rate in the 1.2-19.9 keV range. Indicated are the branches 
in which the source was during the observations. During 
two observation periods (1984 day 251/250, [c], and 1985 
day 232/233, [d]) bursts can be seen which were already re- 
ported by Sztajno et al. (1986). Start times for the different 
observation periods can be found in Table 1. 

Figure 6: 

EXOSAT ME light curve (a) and soft colour curve (b) of 
the observation period 1986 day 093/094. Intensity is the 
dead-time and background corrected count rate in the 1.2- 
19.9 keV range, while the soft colour is the ratio of the count 
rate in the 4.7-6.6 keV and the count rate in the 1.2-4.7 keV 
band. The intensities are not corrected for collimator re- 
sponse (see text). Two new bursts can be seen at (~36 000 
and ~ 104 000 s after the start of the observation). The start 
time of the observation period can be found in Table 1. 

Figure 7: 

High time resolution (0.5 s) argon light curves of the bursts 
which occurred during 1986 day 093 (a) and 1986 day 094 
(b), together with their colour (ratio of the count rates in 
the 4.7-19.9 keV and 1.2-4.7 keV bands) curves, [c] and [d], 
respectively. During the 1986 day 094 burst the observations 
were interrupted for mode changes. The qualified event rate 
from the housekeeping data is therefore also plotted in b 
(with symbol Q) to show the overall burst profile. 

Figure 8: 

(a) A neutron star (radius 10 km) with a symmetric dipole 
field. The spherical magnetosphere (indicated by a dotted 
line) in this plot is taken to be at a radius of 30 km. The last 
closed field line is indicated with a dashed line. The corre- 
sponding accretion polar cap on the surface of the neutron 
star is indicated with a fat line, (b) Same neutron star with 
a shifted (by 5 km) magnetic dipole field. The upper pole 
now has a much larger area than the lower pole. 

Figure Al: 

Observed Poisson levels for the power spectra of HER7 data 
for several LMXBs as a function of the raw observed count 
rate. The HTR4 data of GX 17-1-2 have also been plotted. 
The observed Poisson level was determined by adding it as 
an extra component in the power spectral fits. The different 
symbols refer to the different sources as indicated. Error bars 
are on the order of the symbol sizes. The predicted Poisson 
levels (straight line) and the quadratic fit (dashed line) are 
included in the plot. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science M^^X style file. 
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Table 1: EXOSAT Argon observation log of GX17+2 



Year 


Date 


Day" 


UT start 
(hr:min) 


UT end 
(hr:min) 


Mode'' 


Time 
Resolution 


Spectral 
State 


Config.= 


1983 


Aug 3/4 


215/216 


15:58 


00:36 


E4/T3 


0.3s/8ms 


NB/FB 


H2 


1984 


Sep 5/6 


249/250 


21:47 


05:21 


E4/5/T3'' 


10/0.6s/8ms 


HB 


H1/H2'' 


1984 


Sep 6/7 


250/251 


19:42 


04:04 


E5/T3 


0.6s/8ms 


NB 


Hl/H2= 


1985 


Aug 20/21 


232/233 


10:32 


07:19 


E5/T3 


1 s/ 2 ms 


NB 


H1/H2/WA/ 


1985 


Sep 15/16 


258/259 


16:15 


04:48 


E5/7/T3/49 


0.7s/2/4/0.25ms 


HB 


WA 


1986 


Apr 3/4 


093/094 


04:09 


11:10 


E2/5/7/T4/5'^ 


l,4/ls/8/0.25/lms 


NB/FB 


WA 



"jan 1 = day 1. 

= HER mode = High Energy Resolution mode, T = HTR mode = Higli Time Resolution mode. 
^Array configuration on source: Hl = 'half 1', H2 = 'half 2', WA = whole array. 
'^From 21:47 to 22:07 HER4 (HI) the rest HER5/HTR3 (H2). 
^From 19:41 to 19:51 HI; the rest H2. 

^WA: 10:32 - 16:18, 19:02 - 00:03, 01:52 - 07:19; H2: 16:21 - 17:53, 00:01 - 01:04; HI: 17:57 - 19:02, 01:06 - 01:51. 
9From 16:15 to 03:13 HER7/HTR3 (HER7 effectively one energy band: 1.4-4.1 keV; HTR3: from 16:15 to 22:58 xenon; the 

rest argon); the rest HER5/HTR4. 
''HER2: day 093 15:28-20:31 (Is) and day 093 22:55-day 094 09:07 (4 s); HER5: day 094 09:11-11:10; 

HER7: day 093 04:09-15.14 and day 093 20:36-day 094 09:07; HTR4: day 093 15:25-20:31; 

HTR5: day 093 04:09-15:14 and day 093 20:39-day 094 11:10. 



Table 4: HBO harmonic parameters' 



harmonic 



Obs. 
period 


Sz 




rms 
(%) 


err^ 


Au 
(Hz) 


err 


(Hz) 


err 


rmSH c 
rmSQ 


err 


Ai/H d 


err 


512s 


FFTs 




























1985 


0.116 





.036 


2 


.1 


±0.5 


8 


±3 


49.2 


±0.5 


1.1 


±0.3 


3.0 


±1.4 


day 258 


0.251 





.072 


1 


.3 


±0.4 


4 


±2 


52.1 


±1.0 


0.8 


±0.3 


1.4 


±0.9 


1985 


1.440 





.063 


1 


.2 




5 




60 












day 232 


1.612 





.027 





.5 




5 




60 














1.693 





.018 





.6 




5 




60 












1986 


1.558 





.099 


1 


.1 




5 




60 












day 093 






























16s 


FFTs 




























1985 


1.386 





.027 


1 


.6 




5 




60 












day 232 


1.486 





.045 





.6 




5 




60 














1.603 





.027 


1 


.0 




5 




60 














1.666 





.009 


1 


.0 




5 




60 














1.711 





.009 





.8 




5 




60 












1985 


0.097 





.036 


3 


.0 


±0.4 


16 


±6 


49 


±2 


1.1 


±0.2 


2.9 


±1.3 


day 258 


0.170 





.009 


1 


.4 




6 




52.6 




0.9 










0.287 





.063 


1 


.3 




6 




54.6 




0.9 








HTR4='-'' 


0.079 





.090 


2 


.5 


±0.6 


10 


±8 


52 


±2 


0.7 


±0.2 


1.5 


±1.3 


1986 


1.495 





.090 


1 


.3 




5 




60 












day 093 


1.639 





.027 


1 


.2 




5 




60 














1.738 





.027 





.6 




5 




60 













Whenever a parameter was fixed or not evaluated it is indicated by " — " or "- -" , respectively. 
^AU errors in the power spectral parameters are determined from an error scan through the space 
using Ax^ = 1. If no error is given the parameter value 
is a 1(7 upper limit. 

"^Ratio of the harmonic rms (rmsn) and the HBO rms (rmsg). 
^Ratio of the harmonic FWHM (Ai^n) and the HBO FWHM (Aiyq). 
'^The Poisson level was a free parameter in the fits to these data. 
-'^Effective energy range: 1.4— ll.SkeV. 
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Table 2: Results of power spectral fits of 5f2s FFTs' 



VLFN LFN HFN QFo'' 



Oba. 


Sz 




rms 




ay 


err 


rms 


err 




err 




err 


rms 


err 




err 


rms 


err 




err 


1/ 


err 


'^red 


dof 


period 






% 








% 








(Hz) 




(%) 




(Hz) 




(%) 




(Hz) 




(Hz) 






1983 


1.351 


0.087 


0.63 


±0.16 


0.88 


±0.24 


1.5 


±0.3 


-0.9 


±0.7 


1.2 


±0.9 


1.0 




77 


— 


0.7 




3 


— 


7 


— 


1.0 


36 


day 215 
































— 








— 




— 








1.832 


0.245 


1.01 


±0.08 


1.75 


±0.17 














1.7 


±0.3 


13 


±6 


1.9 




3 


— 


7 


— 


1.2 


39 


1984 


0.584 


0.144 


0.92 


±0.13 


0.79 


±0.11 


3.1 


±0.2 


-1.2 


±0.3 


1.6 


±0.3 


4.5 


±0.4 


147 


±580 


1.2 




3 


— 
— 


31.3 


— 
— 


1.2 


34 


day 249 
































-79 








— 




— 








0.810 


0.054 


0.62 


±0.07 


1.35 


±0.20 


3.1 


±0.2 


-0.7 


±0.2 


2.7 


±0.6 


4.2 


±0.4 


326 


± CO 


1.1 




3 


— 


34.2 


— 


1.0 


34 


































— 250 


















1984 


1.090 


0.072 


0.42 


±0.09 


1.3 


±0.3 


1.4 




— 0.5 


— 


2 


— 


3.7 


±0.3 


38 


±12 


0.7 




3 


— 


7 


— 


1.1 


37 


day 250 




















— 




— 
















— 




— 








1.253 


0.036 


0.72 


±0.09 


0.93 


±0.12 


0.7 




— 0.5 


— 


2 


— 


3.1 


±0.3 


48 


±39 


1.0 




3 


— 


7 


— 


0.8 


37 


1985 


1.440 


0.063 


1.16 


±0.11 


2.17 


±0.21 














2.9 


±0.2 


31 


— 16 
±9 


1.1 




3 


— 


7 


— 


1.4 


40 


day 232 








































— 




— 








1.612 


0.027 


1.13 


±0.05 


1.51 


±0.09 














2.1 


±0.3 


56 


±35 


1.2 


±0.2 


3.0 


±1.6 


7.8 


±0.5 


0.8 


57 




1.693 


0.018 


1.46 


±0.07 


1.65 


±0.09 














2.9 


±0.2 


321 


-18 

± CO 


1.7 


±0.3 


5.2 


±2.4 


6.9 


±0.5 


1.1 


57 


































-178 




















1.756 


0.018 


1.31 


±0.07 


1.62 


±0.11 














3.2 


±0.2 


112 


±82 


1.9 


±0.2 


4.4 


±1.2 


8.1 


±0.3 


1.0 


57 




1.810 


0.018 


2.25 


±0.12 


1.71 


±0.08 














3.1 


±0.2 


269 


— 36 

± CO 


2.1 


±0.1 


3.2 


±0.5 


7.2 


±0.2 


1.5 


57 


































244 




















1.883 


0.036 


1.22 


±0.06 


1.33 


±0.12 














2.4 


±0.6 


108 


± CO 


2.7 


±0.5 


7.3 


±4.4 


6.8 


±0.5 


1.1 


75 


































-91 




-1.3 














1985 


0.116 


0.036 


0.47 


±0.04 


2.2 


±0.5 


3.2 


±0.1 


-0.8 


±0.1 


1.3 


±0.2 


5.4 


±0.3 


20 


±1 


1.9 


±0.1 


2.7 


±0.8 


25.7 


±0.3 


1.0 


105 


day 258 
































-4 




















0.251 


0.072 


0.58 


±0.06 


1.59 


±0.23 


3.6 


±0.1 


-0.6 


±0.1 


1.7 


±0.1 


5.6 


±0.2 


34 


±1 

— 3 


1.6 


±0.2 


2.8 


±0.5 


26.9 


±0.9 


0.8 


102 




0.007 


0.090 


0.53 


±0.08 


1.5 




1.6 


±0.3 


— 1.2 


±0.2 


0.7 


±0.6 


3.5 


±0.5 


18 


±12 


1.3 




3 




24.1 




0.7 


28 






















-1.3 




-0.3 








-6 




















0.296 


0.090 


0.53 


±0.09 


1.42 


±0.20 


1.0 


±0.2 


-3.1 


±1.6 


0.3 


±0.3 


3.1 


±1.8 


30 


±41 


1.7 




3 




27.7 




0.8 


36 






















— 3.3 




— 0.1 




— 0.8 




— 13 


















1986 


1.558 


0.099 


6.9 


±0.5 


2.26 


±0.06 














3.1 


±0.2 


108 


±66 


1.2 




3 




7 




1.4 


40 


day O93S 




















































1.801 


0.072 


6.6 


±0.3 


2.09 


±0.04 














3.2 


±0.2 


41 


±9 


1.3 




3 




7 




1.0 


40 




2.650 


0.298 


5.6 


±0.3 


2.07 


±0.05 














3.0 


±0.2 


40 


±12 


1.2 




7 




16 




1.0 


40 




3.868 


0.406 


6.6 


±0.4 


1.86 


±0.05 














2.4 


±0.3 


127 


-|- CO 

-68 


1.1 




7 




16 




1.3 


40 


HTR4'*>^ 


2.812 


0.162 


6.2 


±0.5 


1.99 


±0.07 














1.7 




77 




1.8 




7 




16 




0.9 


41 



'^Whenever a component was not evaluated it is indicated by -"; whenever we determined an upper limit (icr) of a component, we indicate the errors of the fixed parameters by " — VLFN, LFN and HFN 
are determined from 0.001-1 Hz, 0.01-100 Hz and 0.01-100 Hz, respectively. 
^HBO or NBO. 

"^All errors are determined from an error scan through the space using Ax^ = 1. 
"^The Poisson level was a free parameter in the fits to these data. 
^Effective energy range: 1.4— 4.1keV. 
-^Effective energy range: 1.2— 19.8 keV. 

^Pointing unstable; VLFN results may be corrupted, see text. 



Table 3a: Results of power spectral fits of f6s FFTs (HB/NB) 



VLFN LFN HFN QFo'' 



Obs. 


Sz 




rms 


err ay err 


rms 


err^ 




err 




err 


rms 


err 




err 


rms 


err 


Av 


err 




err 


"red 


dof 


period 






% 




% 








(Hz) 




(%) 




(Hz) 




(%) 




(Hz) 




(Hz) 






1983 
day 215 


1.351 


0.087 


0.5 


1.5 — 


1.8 


±0.2 


-0.2 


±0.3 


2.5 


±1.3 


1.4 




77 


— 


0.9 




3 


4-1 fi 
It I .b 


7 


— 


0.9 


37 


1984 
day 249 


0.495 
0.729 


0.144 
0.027 


1.2 
0.8 


1.5 — 
1.5 — 


3.7 
3.1 


±0.2 
±0.2 


— 0.8 

— 0.8 


±0.3 
±0.2 


1.9 
2.4 


±0.5 
±0.5 


4.6 
4.6 


±0.5 
±0.4 


99 
133 


- 

-57 

-75 

+ CO 

-39 
+ 35 
-15 
±19 


1.2 
1.3 




3 
3 


- 


30.2 
33.2 


- 


1.4 
0.7 


27 
27 




0.838 


0.045 


0.3 


1.5 — 


2.5 


±0.4 


-0.8 


±0.4 


2.0 


±1.0 


4.1 


±0.4 


77 


1.0 




3 




34.5 




1.1 


27 


1984 
day 250 


1.036 


0.054 


0.7 
0.8 


1.5 — 
1.5 — 


1.6 




-0.5 
-0.5 




2 


— 


3.6 


±0.4 
±0.4 


39 
43 


0.4 




3 
3 




7 


— 


1.4 


34 




1.226 


0.018 


0.6 


1.5 — 


0.6 




-0.5 




2 




3.4 


±0.4 


68 


+ 170 

-28 

±6 


0.8 




3 




7 




1.0 


34 




1.280 


0.027 


1.0 


1.5 — 


0.9 




-0.5 




2 




2.6 


±0.4 


19 


1.2 




3 




7 




0.6 


33 


1985 
day 258 


0.097 
0.170 


0.036 
0.009 


1.0 
0.7 


1.5 — 
1.5 — 


2.5 
3.8 


±0.3 
±0.2 


-1.7 
-0.3 


±0.6 
±0.1 


0.7 
2.2 


±0.2 
±0.3 


5.0 
5.8 


±0.2 
±0.2 


10 

32 


±1 
±3 


2.7 
1.5 


±0.3 
±0.2 


5.6 
2.0 


±1.3 
±1.0 


25.1 
26.3 


±0.4 
±0.4 


1.3 
1.0 


44 
88 




0.287 


0.063 


0.9 


1.5 — 


3.5 


±0.1 


-0.6 


±0.1 


1.8 


±0.3 


5.8 


±0.2 


38 


±4 


1.4 


±0.3 


1.7 


±0.9 


27.3 


±0.6 


0.9 


72 


HTR4'*'^ 


0.079 


0.090 


0.4 


1.5 — 


4.8 


±0.1 


-0.1 


±0.1 


3.7 


±0.6 


6.7 




77 




3.4 


±0.3 


6.7 


±2.1 


24.0 


±0.5 


1.2 


42 




-0.0209 


0.072 


1.4 


1.5 — 


1.4 


±0.3 


-0.9 


±1.1 


0.5 


+ 0.9 
-0.3 
±0.1 


3.7 


±0.5 


14 


±5 


1.7 




3 




23.8 




1.1 


30 




0.161 


0.063 


1.0 


1.5 — 


1.4 


±0.2 


-5 


±3 


0.2 


3.1 


±0.5 


13 


±2 


2.0 




3 




26.1 




1.0 


29 




0.359 


0.054 


0.9 


1.5 — 


1.8 




-0.5 




2 




3.6 


+ 1.6 
-0.9 


23 


+ 35 
-13 


1.6 




3 




28.5 




0.8 


33 


^Whenever 


a component 


was not 


evaluated it is indicated by 




wlienever 


we determined 


an upper 


limit (1 


cr) of a 


component, we 


indicate 


tlie errors 


of tlie fixed parameters 


by «— » 


. VLFN, 


LFN and 





HFN are determined from 0.01-1 Hz, 0.01-100 Hz and 0.01-100 Hz, respectively. 
^HBO or NBO. 

""All errors are determined from an error scan through the space using Ax^ = 1- 
"^The Poisson level was a free parameter in the fits to these data. 
'^Effective energy range: 1.4-11.5 keV. 
•^Effective energy range: 1.4-4.1 keV. 

^Negative due to scaling; this indicates that the HB is somewhat longer than the NB. 



Table 3b: Results of power spectral fits of 16 s FFTs (NB/FB) 



Obs. rms err Qy err rms err i^jj err rms err Ai> err v 

period % (%) (Hz) (%) (Hz) (Hz) 



1985 
day 232 



1986 
day 093 



1.386 
1.486 
1.603 
1.666 
1.711 
1.766 
1.801 
1.847 
1.901 
1.495 
1.639 
1.738 
1.856 
2.361 
2.632 
2.993 
3.444 
4.157 
2.659 
3.002 
3.823 



0.027 
0.045 
0.027 
0.009 
0.009 
0.009 
0.009 
0.009 
0.027 
0.090 
0.027 
0.027 
0.045 
0.072 
0.108 
0.099 
0.180 
0.280 
0.126 
0.126 
0.316 



0.7 
0.7 
1.0 
1.0 
0.6 
0.9 
0.8 
0.9 
1.0 
1.0 
0.8 
0.8 
0.9 
1.1 
1.0 
2.6 
2.4 
3.1 
2.1 
2.7 
3.5 



-t- CO 
-0.3 
-1-2.8 
-0.2 



±0.1 
±0.1 



±0.1 

±0.2 

±2.7 
-0.1 
±0.1 

±0.1 

±0.8 
-0.2 
±0.1 

±1.8 
-0.6 
±0.2 

±0.2 

±2.6 
-0.5 
±1.1 
-0.5 
±0.3 



1.5 
1.5 
2.1 
1.8 
1.5 
1.0 
0.9 
1.5 
1.0 
1.6 
1.5 
1.0 
1.3 
1.6 
1.4 
2.5 
1.8 
1.9 
2.2 
2.3 
2.0 



± CO 
— 1.0 
±2.2 
-0.7 



±0.2 
±0.3 



±0.3 

±0.5 

±2.5 
-0.5 
±0.2 

±0.5 

±0.7 

±0.3 

±0.6 

±0.2 

±0.2 

±1.6 
-0.6 
±0.4 

±0.2 



3.3 
2.9 
2.4 
2.8 
2.9 
2.9 
3.3 
3.7 
2.8 
3.0 
3.1 
3.2 
3.4 
3.2 
3.1 
3.5 
2.4 
2.4 
2.5 
1.3 
2.1 



±0.2 
±0.3 
±0.3 
±0.3 
±0.5 
±0.3 
±0.3 
±0.4 
±0.3 
±0.2 
±0.2 
±0.2 
±0.3 
±0.5 
±0.3 
±0.3 
±0.4 
±0.3 
±0.5 
±0.7 



94 
108 
63 
82 
70 
160 
141 
33 
62 
82 
109 
149 
105 
131 
62 
52 
151 
132 
12 



±36 
±43 

±81 

— 36 
±42 

±94 

-31 

±116 

-49 

±73 

-36 

±29 

— 12 
±77 
-28 
±21 

±29 

±34 

±33 

±54 

±19 

±18 

±165 

-74 

±64 



± CO 
-11 



1.3 
1.5 
1.3 
1.2 



1.9 
2.5 
1.2 
1.0 
0.9 
1.5 
2.3 
1.3 
1 .2 
1.1 
0.5 
1.0 



±0.3 



±0.8 
±0.2 
±0.2 



±0.3 
±0.3 



3 
3.9 

3 

3 

6 
3.4 
2.3 

3 
4.1 

3 

3 

3 
2.7 

7 

7 

7 

7 

7 

7 



±2.0 



±6 

±1.2 

±0.7 



±1.1 
±3 



7 
5.3 
7 
7 
6.8 
8.1 
7.6 
7 
7.4 
7 
7 
7 
6.9 
16.1 
16 
16 
16 
16 
16 



±0.7 



±0.9 
±0.4 
±0.2 



±0.3 
±0.6 



1.3 
1.0 
0.8 
1.1 
1.0 
1.5 
0.8 
1.8 
0.9 
1.0 
1.2 
1.4 
0.8 
1.0 
0.8 
1 .0 
0.9 
0.8 
1.2 
1.3 
1.5 



30 

53 



28 

53 



31 
31 
31 
57 
57 
31 
31 
31 
31 
24 



^Whenever a component was not evaluated it is indicated by -"; wlienever we determined an upper limit (Icr) ol a component, we indicate tlie errors ol 

the fixed parameters by " — ". VliFN and HFN are determined from 0.01-1 Hz and 0.01-100 Hz, respectively. 
''NBO or FF)0. 

^All errors are determined from an error scan through the space using A^*^ — 1. 

^The Poissou level was a free parameter in the fits to these data. Effective energy range: 1.2— 19.8 keV. 



5 



Table 5: Upper limits to 
modulation strength 



Burst 


(Hz) 


Upper 
limit (% )*> 


I 


64 


6.1 


II 


256 


1.0 


III 


512 


2.1 


IV 


512 


1.9 


II+III+IV 


256 


0.9 


III+IV 


512 


1.5 



^i^Ny is the Nyquist frequency. 
^99 per cent confidence upper limit 

to the fractional sinusoidal amplitude 

of the total flux. 



Table 6: Polar cap areas (km^) 



Dipole field offset (km) 



(km) 





1 


2 


3 


4 


5 


15 


259 


325 


389 


456 


519 


591 




259 


214 


157 


109 


70 


41 


40 


84 


110 


140 


172 


207 


245 




84 


68 


49 


33 


22 


13 


100 


33 


43 


55 


69 


86 


106 




33 


26 


19 


13 


8 


5 



"^Total area of the neutron star is ~1257km^. 
^Trn is the magnetospheric radius. 



Table Al: EXOSAT Argon HER7 observation log 



Source 


Year 


Date 


Day" 


Time 
Resolution 


Spectral 
State 


GX5-1 


1985 


Aug 27/28 


239/240 


3.9ms 


NB/FB 




1985 


Aug 29 


241 


3.9ms 


NB 




1985 


Aug 30 


242 


3.9ms 


NB 




1985 


Sep 17 


260 


3.9ms 


HB 


CygX-2 


1985 


Oct 28/29 


301/302 


3.9ms 


HB 




1985 


Oct 29 


302 


3.9ms 


HB 




1985 


Nov 14/15 


318/319 


3.9ms 


NB/FB 


GX17+2 


1986 


Apr 3/4 


093/094 


7.8 ms 


NB/FB 


ScoX-1 


1985 


Aug 24 


236 


3.9ms 


FB 




1985 


Aug 25 


237 


3.9ms 


NB/FB 




1986 


Mar 11/12 


070/071 


7.8 ms 


HB/NB 




1986 


Mar 13/14 


072/073 


3. 9/7. 8ms 


NB 


GX3+1 


1985 


Sep 4 


247 


3.9ms 


UB*" 


^Jan 1 = day 1. 
^UB = upper banana 


(see Hasinger & 


van der Klis 


1989). 
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